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ABSTRACT 


This  newly  initiated  project  will  develop  and  test  a  methodology  for  locating  seismic  events  from  combined  data  sets 
of  regional  and  teleseismic  arrival  times,  based  on  consistent  travel-time  predictions  from  a  unified  3D  Earth  model. 
One  focus  of  the  project  is  to  address  the  practical  difficulty  of  raytracing  in  3D  models,  which  has  been  a  serious 
impediment  to  the  pursuit  of  3D  event  location  methods.  We  will  investigate  whether,  for  teleseismic  travel-time 
prediction,  approximate  techniques,  in  particular  linearization  around  rays  calculated  in  a  ID  reference  model,  are 
adequate  for  the  purpose  of  event  location  as  they  are  commonly  assumed  to  be  for  the  purpose  of  global 
tomography.  A  second  focus  of  the  project  is  tomographic  calibration  of  a  3D  model  with  combined  regional  and 
teleseismic  data  from  earthquakes  and  ground-truth  events.  We  will  investigate  the  hypothesis  that  a  joint 
regional/teleseismic  calibration  will  lead  to  a  noticeable  improvement  in  location  accuracy  over  the  modest  and 
inconsistent  improvements  3D  models  have  yielded  to  date.  Additionally,  we  will  consider  whether  travel-time 
prediction  errors  inferred  from  a  tomographic  uncertainty  analysis  can  provide  an  appropriate  weighting  of  various 
teleseismic  and  regional  phases  to  optimize  location  accuracy  even  further.  We  plan  to  validate  our  methodology 
with  catalog  data  from  Asia,  the  Middle  East  and  North  Africa,  using  available  regional  crust/upper  mantle  models 
for  these  areas  in  conjunction  with  published  3D  global  models  of  the  deeper  mantle.  Special  attention  will  be  given 
to  south-central  Asia,  where  Weston  Geophysical  and  MIT  are  currently  applying  body-wave  and  surface-wave 
tomography  to  develop  a  regional  model  of  the  crust  and  upper  mantle.  This  project  will  attempt  to  improve  and 
extend  the  Weston/MIT  model  by  adding  teleseismic  constraints. 
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OBJECTIVE 


The  development  of  three-dimensional  models  of  Earth  strueture  continues  to  be  an  active  field  of  seismology 
advancing  along  several  fronts.  A  small  sampling  of  these  efforts  include  global  travel-time  tomography  (e.g.,  van 
der  Hilst  et  al.,  1997;  Antolik  et  al.,  2003),  global  surface-wave  dispersion  tomography  (e.g.,  Ritzwoller  et  al., 

2002),  regional  models  for  Asia  (Murphy  et  al.,  2005;  Li  et  al.,  2006)  and  the  Middle  East  (Pasyanos  and  Walter, 
2002).  An  important  question  for  nuclear  monitoring  is  whether  these,  or  improved,  3D  Earth  models  can  produce 
routine  event  locations  with  significantly  reduced  errors  compared  to  ID  models  when  data  at  all  event-station 
distances  are  used  together. 

The  answer  to  this  question  would  seem  to  be  a  definite  yes  when  you  consider  that  travel-time  predictions  from  a 
ID  model  like  akl35  (Kennett  et  al.,  1995)  differ  from  observed  travel  times  systematically  by  up  to  10  seconds  at 
regional  distances  and  a  few  seconds  at  teleseismic  distances,  more  than  enough  to  induce  event  mislocations  of  tens 
of  kilometers.  An  example  from  our  own  work  on  tomography  and  location  supports  this  view.  Weston  Geophysical 
Corp.  and  MIT  are  developing  a  new  3D  cmst  and  upper  mantle  model  for  southern  Asia  (see  Figure  1)  based  on  Pn 
travel  times  from  the  EFIB  database  (Engdahl  et  al.,  1998)  and  Rayleigh-wave  dispersion  in  the  period  range  10-150 
sec  (see  Reiter  and  Rodi,  2008).  The  inversion  model,  which  we  call  JWM,  constrains  the  P  and  S  velocities  only  at 
depths  less  than  410  km,  so  we  are  validating  JWM  hy  relocating  ground-truth  events  with  local  and  regional 
arrival-time  data.  The  results  of  two  such  tests  are  shown  in  Table  1,  where  mislocations  are  listed  for  the 
1998/05/05  nuclear  test  in  Pakistan  and  the  1998/05/28  nuclear  test  in  India.  We  see  that,  even  with  a  large  number 
of  regional  phases,  the  locations  based  on  akl35  are  in  considerably  more  error  than  the  locations  derived  from  our 
3D  model.  (These  tests  used  only  regional  P  arrivals  (Pn)  and  the  results  shown  are  for  constrained  focal  depth.) 


Figure  1.  Study  regiou  for  the  Westou/MIT  regioual  tomography  project,  curreutly  uuderway.  The  project  is 
applyiug  joiut  hody-wave  travel-time  aud  surface-wave  dispersiou  tomography  to  derive  a  self- 
cousisteut  model  for  the  P-wave  aud  5-wave  velocities  of  the  crust  aud  upper  mautle.  The  iuversiou 
model  {JWM)  is  descrihed  hy  Reiter  aud  Rodi  (2008). 
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Table  1:  Location  errors  for  the  1998  Pakistan  and  India  nnclear  tests,  nsing  regional  arrival-time  data. 


Event 

Model 

No.  of  data 

Epicenter 

(for  reg.) 

Error  (km) 

Pakistan  test 

akl  3  5 

35 

12.6 

Pakistan  test 

JWM 

35 

2.9 

India  test 

akl35 

27 

37.2 

India  test 

JWM 

27 

16.3 

The  results  in  Table  1  are  very  eneouraging  sinee  events  of  monitoring  interest  are  typieally  small  and  yield  mostly 
regional  observations.  However,  sueh  events  often  yield  teleseismie  observations,  too,  as  in  the  ease  of  the  North 
Korea  test  on  2006/10/08.  Furthermore,  the  loeation  of  large  events  is  also  important  for  nuelear  monitoring  sinee 
well-loeated  events,  of  any  size,  serve  as  referenee  events  for  ealibration  and  modeling  studies.  Therefore,  we  also 
reloeated  the  explosions  in  Pakistan  and  India  with  the  addition  of  the  many  teleseismie  P-wave  observations 
available,  using  stations  up  to  50°  from  the  events.  We  still  used  JWM  for  predieting  Pn  travel  times  but,  sinee  our 
model  is  eonstrained  only  to  410  km  in  depth,  we  used  akl35  for  predieting  teleseismie  travel-times.  We  eompared 
this  mixed  approaeh  with  purely  akl35  solutions.  The  results  are  in  Table  2. 

Table  2:  Location  errors  for  the  1998  Pakistan  and  India  nnclear  tests,  nsing  regional  and  teleseismie  data. 


Event 

Model 

No.  of  data 

Reg.  Tel. 

Epicenter 
Error  (km) 

Pakistan  test 

akl  3  5 

35 

85 

13.0 

Pakistan  test 

JWM 

35 

85 

16.6 

India  test 

akl  3  5 

27 

99 

15.2 

India  test 

JWM 

27 

99 

5.9 

When  akl35  is  used  for  both  regional  and  teleseismie  times,  the  epieenter  misloeation  of  the  Pakistan  test  is  about 
the  same  with  or  without  teleseismie  data  (line  1  of  Tables  1  and  2).  However,  when  JWM  is  used  for  regional  times 
and  akl35  for  teleseismie,  the  misloeation  inereases  signifieantly  (line  2  of  the  tables).  This  degradation  does  not 
happen  in  the  ease  of  the  India  explosion,  however  (lines  3  and  4  of  the  tables),  for  whieh  the  addition  of  teleseismie 
data  reduees  the  misloeation  with  either  regional  model.  Notably,  JWM  still  performs  better  than  akl35.  The 
eonelusion  we  draw  is  that  mixing  a  ealibrated  regional  model  with  akl  3 5  (whieh  is  only  ealibrated  in  an  average 
global  sense)  ean  make  event  loeations  worse,  although  not  every  time  as  we  have  seen. 

Other  studies  have  demonstrated  improvement  in  regional  event  loeations  based  on  3D  erust/upper  mantle  models, 
although  they  have  been  restrieted  mostly  to  a  priori,  rather  than  tomographie,  models  (Ryaboy  et  al.,  2001;  Johnson 
and  Vineent,  2002;  Pasyanos  et  al.,  2004;  Flanagan  et  al.,  2007).  However,  relatively  little  work  has  been  on  joint 
regionaFteleseismie  loeation  with  3D  models.  In  one  of  the  few  published  studies,  Yang  et  al.  (2004)  found  that 
event  loeations  obtained  by  eombining  regional  and  teleseismie  data  were  only  more  aeeurate  than  loeations  found 
from  the  separate  data  sets  if  they  eorreeted  for  a  bias  between  the  teleseismie  and  regional  travel-time  predietions, 
eaeh  made  from  a  different  3D  model  (CUBl.O  for  regional,  J362  for  teleseismie).  The  eorreetion  they  inferred  is 
small  (0.79  see),  whieh  is  eonsistent  with  the  faet  that  3D  erust/upper  mantle  effeets  on  teleseismie  travel  times  are 
relatively  small  (<  2  see).  But  event  small  systematie  shifts  in  teleseismie  travel  times  ean  induee  large  loeation 
errors  owing  to  the  high  apparent  veloeity  of  teleseismie  arrivals. 

It  is  elear  to  us  that  the  most  promising  approaeh  to  event  loeation  with  teleseismie  and  regional  observations  is  to 
eomputer  travel  times  eonsistently  from  a  single  model.  Further,  that  model  should  be  jointly  ealibrated  with  both 
regional  and  teleseismie  data. 

RESEARCH  PLANNED 


This  projeet  will  develop  and  test  a  self-eonsistent  unified  methodology  for  seismie  event  loeation  and  loeation 
ealibration  based  on  3D  Earth  models.  Some  of  the  questions  we  will  address  are  the  following: 
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■  Given  a  unified  3D  model  of  the  Earth,  what  numerical  algorithms  are  best  suited  to  the  forward  modeling 
of  travel  times  for  various  seismic  phases  for  use  in  event  location?  The  phases  we  will  be  concerned  with 
are  the  local,  regional  and  teleseismic  phases  that  are  routinely  available  in  earthquake  catalogs  such  as  the 
EHB  database,  and  those  potentially  observable  in  monitoring  applications. 

■  Using  observations  of  the  same  phases  from  reference  events  and  earthquake  catalogs,  or  from  temporary 
station  deployments,  what  is  the  best  strategy  for  calibrating  the  unified  Earth  model?  An  important  issue 
here  is  whether  regional  models  of  the  crust  and  upper  mantle  can  be  calibrated  separately  from  lower 
mantle  models,  or  whether  a  joint  regional/teleseismic  tomography  is  essential.  Related  to  this,  can  a 
tomographic  inversion  for  the  crust  and  upper  mantle  include  constraints  from  teleseismic  data  without 
varying  the  deeper  mantle  structure? 

■  Can  a  rigorous  approach  to  model  uncertainty  provide  automatic  data  weights  for  an  event  location 
algorithm,  such  that  well-calibrated  and  weakly  or  un-calibrated  travel-time  predictions  are  combined 
optimally?  Location  algorithms  used  in  nuclear  monitoring  currently  employ  estimates  of  “model  errors”  to 
perform  this  weighting.  We  will  investigate  whether  more  formal  uncertainty  estimates  derived  from 
tomographic  calibration,  including  covariances  between  travel-time  predictions  along  similar  paths,  can 
perform  better. 

We  plan  to  address  these  questions  using  a  variety  of  available  3D  Earth  models,  including  the  JITM  model  referred 
to  above,  the  WENAl.O  model  developed  at  LENT  (Pasyanos  et  al.,  2004),  and  global  mantle  models  such  as 
J362D28  (Antolik  et  al.,  2003)  and  MITP07  (Li  et  al.,  2008). 

Statement  of  the  Problem 

The  forward  model  in  seismic  event  location  employs  a  set  of  travel-time  functions,  each  predicting  the  travel  time 
of  a  seismic  phase  from  any  potential  event  hypocenter  to  an  observing  station.  This  prediction  is  based  on  wave 
propagation  through  a  reference  velocity  model.  Since  the  Earth's  true  velocity  differs  from  the  reference  model,  the 
forward  model  admits  a  prediction  error  (Billings  et  al.,  1994). 

Let  Vo  denote  the  velocity  of  the  reference  model  (focusing  on  just  the  velocity  for  P  waves),  and  let  Ve  be  the  Earth's 
true  velocity.  We  can  express  the  prediction  error  as 

e(x)  =  T(\;ve)  -  T{\;vo),  (1) 

where  the  travel-time  function,  T{\;v),  integrates  1/v  along  the  raypath  connecting  the  station  to  the  hypocenter  x. 
Prediction  errors  add  to  the  observational  errors  in  arrival-time  picks,  and  a  difficult  but  important  problem  is  to 
account  for  both  kinds  of  error  in  locating  an  event. 

Calibration  is  the  task  of  reducing  or  correcting  the  prediction  errors  in  travel  times.  One  calibration  approach  that 
has  been  used  in  nuclear  monitoring  applications  (e.g.,  Myers  and  Schultz,  2000)  is  the  kriging,  or  geostatistical 
interpolation,  of  observed  travel-time  residuals  from  ground-truth  events,  which  leads  to  a  correction  function  or 
surface,  A7’(x),  for  a  given  station/phase.  This  method  works  well  when  all  of  the  network  stations  have  observed  a 
GT  event  (or  events)  near  an  event  of  interest.  The  forward  model  used  in  event  location  is  then  taken  as  T{\;vo)  + 
A7’(x).  Model-based  calibration,  in  contrast,  attempts  to  correct  the  velocity  model  itself:  replacing  Vo  with  an 
improved  model  Vc  -  Vo  -l-  Av.  The  corrected  model  can  be  derived  with  travel-time  tomography,  applied  to  data 
similar  to  that  used  in  kriging,  but  also  can  incorporate  a  variety  of  other  constraints.  An  advantage  of  model-based 
calibration  is  that  it  allows  the  calibration  of  stations  and  seismic  phases  which  have  not  been  directly  observed,  as 
long  as  Vc  is  defined  for  relevant  portions  of  the  Earth  and,  additionally,  appropriate  ray  tracing  techniques  are 
available  for  the  phases  of  interest. 

The  conceptual  goal  of  this  project  can  now  be  stated  as  follows.  Given  a  globally  defined  3D  model,  Vc,  we  wish  to 
identify  forward  modeling  and  event  location  techniques  which  allow  us  to  take  advantage  of  all  the  travel-time  data 
available  for  locating  an  event.  Conversely,  we  wish  to  develop  tomographic  methods  for  calibrating  Vc,  given  the 
available  travel-time  observations  from  ground-truth  and  other  past  events. 
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3D  Raytracing 

A  critical  choice  in  both  event  location  and  travel-time  tomography  with  3D  Earth  models  is  the  choice  of  numerical 
technique  to  use  for  calculating  travel  times  between  sources  and  receivers.  The  choice  involves  somewhat  different 
issues  in  location  and  tomography.  First,  an  event  locator  requires  travel-time  calculations  for  arbitrary  hypocenters 
as  it  searches  for  a  best-fitting  event  location,  while  tomography  is  usually  performed  with  event  locations  held 
fixed.  Second,  tomography  needs  the  sensitivity  of  each  calculated  travel  time  to  the  velocity  structure  of  the  Earth, 
i.e.  the  ray  trajectory  connecting  source  and  receiver,  or  its  equivalent. 

For  travel-time  calculation  in  a  ID  reference  model  (e.g.,  akl35)  the  choice  is  straightforward  since  analytical 
solutions  for  travel  times  and  rays,  for  any  seismic  phase  (direct,  refracted,  diffracted  or  converted)  are  available  and 
easily  implemented  in  fast  algorithms  (Buland  and  Chapman,  1983).  In  contrast,  travel-time  calculation  in  3D 
models  presents  many  challenges. 

Global  travel-time  tomography  has  generally  been  performed  with  a  linearized  approximation  to  the  travel-time 
function.  In  this  approximation,  the  slowness  of  the  3D  model  (1/vc)  is  integrated  along  the  raypaths  computed  in 
the  reference  model,  Vq,  with  the  reference  model  taken  as  ID.  We  can  identify  the  ray  connecting  a  given  station 
and  a  given  event  hypocenter,  Xev,  with  a  sensitivity  kernel,  a(x,Xev;vo),  that  is  concentrated  along  the  ray  trajectory. 
Then  the  linearized  approximation  to  the  travel  time  through  vc  can  be  written  as 

7’(Xev;vc)  =  7’(Xev;vo)  + 1  dx  a(x,Xev;vo)  [Vc '(x)  -  Vo"'  (x)] 

=  1  dx  a(x,Xev;vo)  Vc '(x).  (2) 

We  point  out  that  this  equation  also  applies  to  finite-frequency  travel  times,  which  Li  et  al.  (2006)  and  others  have 
used  to  model  long-period  body-wave  observations.  In  this  case,  a(x,Xev;vo)  becomes  a  “banana-doughnut” 
sensitivity  kernel  (e.g.,  Dahlen  et  al.,  2000)  evaluated  for  the  reference  velocity  model. 

For  the  calculation  of  local  and  regional  direct  P  arrivals,  linearization  is  not  adequate,  as  we  will  see  in  a  moment. 
For  this  reason,  Reiter  and  Rodi  (2006)  use  the  finite-difference  ray  fracing  method  due  to  Podvin  and  Lecomte 
(1991),  as  implemented  in  software  by  Lomax  (2000).  Reiter  et  al.  (2005)  augmented  the  Podvin-Lecomte  (P-L) 
algorithm  with  a  back-chaining  algorithm  to  generate  sensitivities  of  travel  times  with  respect  to  the  velocity 
parameters  of  the  model.  Since  the  P-L  algorithm  is  written  for  flat-Earth  models  in  Cartesian  coordinates,  we  have 
integrated  it  into  our  tomography  software  with  appropriate  transformations  between  model  parameterizations  and 
sensitivities.  Reiter  et  al.  (2005)  provide  additional  details  of  the  ray  tracing  techniques  we  use  for  regional 
tomography. 

We  illustrate  the  importance  of  3D  ray  tracing  in  modeling  regional  travel  times  using  our  JWM  tomographic  model. 
First,  we  compare  travel  times  calculated  through  JWM,  using  the  Podvin-Lecomte  algorithm,  to  the  travel  times 
calculated  through  akl35.  In  each  case,  the  times  are  calculated  for  the  76,000  event-station  pairs  in  our  tomography 
data  set  (extracted  from  the  EFIB  catalog).  The  difference  between  the  JWM  and  akl35  times  for  these  paths  is 
plotted  as  a  function  of  epicenfral  distance  in  the  left  panel  of  Figure  2.  Not  surprisingly,  the  models  predict  Pn 
times  differing  by  several  seconds. 

Can  Pn  travel  times  through  JWM  he  approximated  with  linearization  around  akl35,  using  Equation  (2)?  The  right 
panel  of  Figure  2  shows  the  difference  between  the  P-L  calculated  times  (“Nonlin.”)  and  the  akl35  linearized 
approximation  (“Linear”).  We  see  that  the  errors  in  the  linearized  times  are  as  high  as  25  sec.  Much  of  this 
difference  owes  to  the  thick  crust  over  much  of  the  model  region,  exceeding  75  km  thickness  beneath  the 
Flimalayas.  Therefore,  large  pieces  of  the  akl35  Pn  raypaths  are  incorrectly  located  within  the  JIFM  crust.  While  it 
might  be  possible  to  correct  for  this  effect,  or  optimize  the  ID  model  used  to  obtain  the  reference  raypaths,  the  fact 
is  that  integrating  a  3D  velocity  model  along  reference  rays  from  a  ID  model  cannot  accurately  model  Pn  travel 
times  through  a  complex,  3D  Earth  model. 
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Fignre  2.  Experiments  with  regional  travel-time  calcnlation  performed  with  the  3D  tomographic  model  JWM 
of  Reiter  and  Rodi  (2006).  Left:  Difference  between  calcnlated  travel  times  for  JWM  and  akl35. 
The  JWM  times  were  calcnlated  with  the  Podvin-Lecomte  (P-L)  finite-difference  algorithm.  Right: 
Difference  between  calcnlated  travel  times  from  JWM,  obtained  two  ways:  “nonlinear”  times 
obtained  by  the  P-L  method,  and  “linearized”  times  obtained  by  integrating  the  JWM  slowness 
over  akl35  raypaths.  In  both  panels,  times  are  plotted  for  the  76,000  Pn  paths  nsed  in  the  JWM 
tomography. 

However,  it  is  interesting  that  the  linearized  travel  times  are  not  too  bad  for  the  longest  Pn  paths,  i.e.,  beyond  14°, 
corresponding  to  rays  bottoming  between  210  and  410  km  in  the  akl35  model.  This  is  because  the  crustal  thickness 
variations  affect  these  raypaths  much  less  than  for  shorter  epicentral  distance,  making  Equation  (2)  more  accurate. 
Also,  the  velocity  differences  between  JWM  and  akl35  in  the  mantle  are  smaller  than  in  the  crust. 

Therefore,  it  is  reasonable  to  ask  whether  the  linearization  approach  would  even  more  accurately  predict  travel  times 
in  the  teleseismic  distance  range.  This  is  a  question  we  will  investigate  in  detail  as  part  of  the  proposed  project.  A 
preliminary  look  into  the  matter  is  offered  in  Figure  3.  This  figure  is  analogous  to  Figure  2,  but  the  travel-time 
calculations  were  performed  only  for  a  line  of  stations  emanating  from  station  JYP  (33°N,  75°E)  along  a  great  circle 
at  an  azimuth  of  -30°  from  north,  out  to  near-teleseismic  distances.  To  calculate  teleseismic  times  in  JWM  with 
Podvin-Lecomte,  JWM  was  extended  below  410  km  depth  with  the  akl35  velocity  profile.  The  left  panel  of  Figure  3 
shows  that  the  influence  of  the  3D  stmcture  in  JWM  on  teleseismic  times  is  less  than  1  second,  given  this 
extrapolation  at  depth.  The  panel  on  the  right,  nonetheless,  shows  that  the  teleseismic  travel  times  through  JWM  are 
accurately  modeled  by  the  linearized  approach  applied  with  akl35  raypaths.  It  remains  to  be  seen  whether  the 
presence  of  deep  3D  structure  alters  this  conclusion. 


Figure  3.  Left:  Difference  between  calculated  travel  times  from  JWM  and  akl35.  Right:  Calculation  error  of 
travel  times  from  JWM  by  using  linearization. 
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Event  Location  with  a  3D  Model 

The  location  validation  effort  in  the  Weston/MIT  tomography  project  for  southern  Asia  is  using  the  grid-search 
location  algorithm  known  as  GMEL  (Grid-search  Multiple-Event  Location),  developed  at  MIT  under  previous 
projects  and  extended  to  3D  models  in  a  previous  joint  project  with  Weston  Geophysical  Corp.  GMEL  uses 
travel-time  predictions  for  either  ID  reference  velocity  models  or  3D  reference  models,  using  interpolation  of 
pre-calculated  tables  in  each  case.  A  summary  description  of  the  numerical  techniques  used  in  GMEL  is  available  in 
Rodi  (2006).  We  plan  to  use  GMEL  in  the  proposed  work,  but  we  will  investigate  whether  teleseismic  travel  times 
for  ID  models  and  linearized  travel  times  for  3D  models  can  be  efficiently  calculated  during  the  grid-search  location 
process  (“on-the-fly”).  This  will  address  an  important  implementation  issue  in  the  use  of  3D  models  for  event 
location  and,  if  on-the-fly  calculations  are  feasible,  GMEL  will  be  a  more  useful  tool  for  the  project. 

Another  previous  effort  we  plan  to  take  advantage  of  in  the  proposed  work  is  a  concurrent  project  on  covariance 
modeling.  Rodi  and  Myers  (2006)  presented  a  method  for  computing  variances  of,  and  covariances  between, 
travel-time  predictions  from  a  reference  model,  based  on  a  geostatistical  description  of  the  velocity  anomalies  in  the 
Earth  relative  to  the  reference  model.  The  numerical  technique  for  doing  this  integrates  travel-time  sensitivities  with 
a  velocity  model  covariance  function,  C(x,x'),  that  characterizes  the  velocity  heterogeneity.  Thus,  the  covariance 
between  the  travel-time  predictions  along  two  paths,  i  and  j,  is  given  by 

(T,y(Xev)  ^\dx\  d\'  a,(x,Xev;vo)  Uj  (x,Xev;vo)  C(x,x'),  (3) 

where  a,  and  a,  are  the  sensitivity  kernels  for  the  two  paths.  An  example  of  this  calculation,  taken  from  Rodi  and 
Myers  (2006),  is  shown  in  Figure  4. 

We  plan  to  implement  the  use  of  travel-time  covariances  in  GMEL,  and  to  extend  the  technique  by  replacing  C(x,x') 
with  the  a  posteriori  covariance  that  results  from  performing  travel-time  tomography.  This  capability  will 
accomplish  a  key  part  of  our  unified  approach,  with  the  data  in  an  event  location  automatically  and  optimally 
weighted  in  correspondence  to  how  well  their  travel-time  predictions  are  calibrated. 

CONCLUSIONS  AND  RECOMMENDATIONS 


Replacing  a  ID  Earth  model  with  a  global  3D  Earth  model  for  routine  event  location  holds  the  promise  of  achieving 
much  greater  location  accuracy,  but  it  is  an  enormous  undertaking  when  one  considers  all  of  the  steps  involved: 
model  integration,  calibration  and  validation,  and  efficient  implementation  for  use  in  a  locator.  Moreover,  it  is  an 
undertaking  that  is  highly  vulnerable  to  the  one-step-forward,  two-steps-back  syndrome,  with  improvement  on  one 
front  (e.g.,  better  regional  travel-time  predictions)  offset  by  regression  on  another  (e.g.,  discrepancy  between 
regional  models,  inconsistency  with  teleseismic  predictions).  The  greatest  benefit  of  our  project  will  be  a  strategy  for 
achieving  the  maximum  effect  from  future  efforts  to  develop  and  implement  3D  Earth  models  for  event  location. 
Additionally  we  expect  this  project  to  have  some  concrete  benefits  beyond  the  methodology  it  develops.  These 
include  an  improvement  to  the  Weston/MIT  model  of  the  crust  and  upper  mantle  in  southern  Asia,  incorporating 
teleseismic  constraints,  and  an  improved  version  of  the  EHB  catalog  for  southern  Asia. 
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Fignre  4.  Top  -  Diagram  of  the  network  geometry  nsed  for  testing  onr  nnmerical  approach  to  modeling 

travel-time  covariances  among  a  set  of  stations.  There  are  12  stations  at  regional  distances  from 
an  earthqnake  in  northeastern  Iran.  Bottom  Left  -  The  standard  deviation  of  the  travel-time 
prediction  error  at  each  of  the  12  stations,  plotted  as  a  fnnction  of  epicentral  distance.  Bottom 
Right  -  Correlation  coefficient  between  station  pairs,  plotted  versns  inter-station  azimnth.  The 
calcnlations  assnmed  that  velocity  variations  in  the  Earth  have  a  standard  deviation  of  2%  and 
correlation  lengths  of  300  km  laterally  and  100  km  vertically. 
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